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The nature of Calcium-Silicate-Hydrate (C-S-H), the binding phase of cement, remains a con¬ 
troversial question. In particular, contrary to the former crystalline model, it was recently proposed 
that its nanoscale structure was actually amorphous. To elucidate this issue, we analyzed the struc¬ 
ture of a realistic simulation of C-S-H, and compared the latter to crystalline tobermorite, a natural 
analogue to cement, and to an artificial ideal glass. Results clearly support that C-S-H is amor¬ 
phous. However, its structure shows an intermediate degree of order, retaining some characteristics 
of the crystal while acquiring an overall glass-like disorder. Thanks to a detailed quantification of 
order and disorder, we show that its amorphous state mainly arises from its hydration. 


INTRODUCTION 

Calcium-silicate-hydrate (C-S-H) is the main hydra¬ 
tion product in Portland cement and acts as a binding 
phase in the paste |T| . Thus, it is responsible for its 
strength, durability and creep properties [2]. 

Despite its ubiquitous presence in the built environ¬ 
ment, the structure of C-S-H at the nanoscale remains 
controversial. As X-Ray diffraction patterns from C-S- 
H have shown to exhibit only a few broad and weak 
diffraction maxima, it has been described as an amor¬ 
phous material mg. However, most experimental stud¬ 
ies mm suggest that its structure is close to the one of 
tobermorite, although their composition differ. Hence, 
it is still unclear whether C-S-H should be considered 
as a crystalline or an amorphous material. Fortunately, 
a realistic atomistic model of C-S-H has recently been 
reported [ 11 1 - ITJ] . thus opening the way to elucidate its 
nature. 

To quantify the atomic order in C-S-H, we analyzed 
the previously mentioned atomistic model. We compared 
its structure with the one of tobermorite, its crystal ana¬ 
logue, while taking care of rescaling the results to account 
for the difference of chemistry. The structure was also 
compared to the one of an ideal artificial glass, formed by 
quickly heating and cooling a C-S-H system. By doing 
a consistent analysis of the structure of those three sys¬ 
tems, simulated in the same conditions and with the same 
potential, we were able to assess their relative atomic or¬ 
der both at the short and medium range order. 

The article is organized as follows. In section [] we 
present the numerical model and the methodology that 
we used to simulate C-S-H and the corresponding crystal 
and glass. Then, we analyze the atomic structure of C- 


S-H as compared to the ones of the crystal and the glass 
in the short range oder in [] This comparison is then 
extended to the medium range order in [] Concluding 
remarks are presented in secion [] 

METHODS 

In this section, we detail the model and the procedure 
used to simulate C-S-H and its equivalent crystal and 
glass. 

A realistic model of C—S—H 

To describe the disordered molecular structure of C- 
S-H, Pellenq et al. jT5] proposed a realistic model for C- 
S-H with the stoichiometry of (CaO)i . 65 (Si0 2 )(H 2 0) 1 . 73 - 
We generated the C-S-H model by introducing defects in 
an 111 tobermorite m configuration, following a combi¬ 
natorial procedure. llA tobermorite consists of pseudo- 
octahedral calcium oxide sheets, which are surrounded by 
silicate tetrahedral chains. The latter consists of bridging 
oxygen atoms and Q 2 silicon atoms (having two bridg¬ 
ing and two non-bridging terminal oxygen atoms) [ T7] . 
Those negatively charged calcium-silicate sheets are sep¬ 
arated from each other by an interlayer spacing, which 
contains interlayer water molecules and charge-balancing 
calcium cations. While the Ca/Si ratio in llA tober¬ 
morite is 1, this ratio is increased to 1.65 in the present 
C-S-H model through randomly removing Si0 2 groups. 
The defects in silicate chains provide possible sites for 
adsorption of extra water molecules. The adsorption of 
water molecules in the structurally defected tobermorite 
model was performed via the Grand Canonical Monte- 
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Carlo method, ensuring equilibrium with bulk water at 
constant volume and room temperature. The REAXFF 
potential j!8| . a reactive potential, was then used to ac¬ 
count for the reaction of the interlayer water with the 
defective calcium-silicate sheets. The use of the reac¬ 
tive potential allows observing the dissociation of wa¬ 
ter molecules into hydroxyl groups. More details on the 
preparation of the model and on the multiple validations 
with respect to experiments can be found elsewhere m 

m. 


Simulation of C—S H 

We simulated the C-S-H model previously presented, 
made of 501 atoms, by molecular dynamics using the 
LAMMPS package m■ To this end, we used the 
REAXFF potential jOE] with a time step of 0.25fs. We 
first relaxed the system at zero pressure and 300K dur¬ 
ing 2.5ns in the NPT ensemble and made sure that the 
convergence of the energy and volume was achieved. We 
then run a 25ps simulation in the NVT ensemble for sta¬ 
tistical averaging. 


Simulation of the crystal analogue 

We chose to study llA tobermorite as it is considered 
as a natural crystal analogue for C-S-H (see section []). 
This choice was also motivated by the fact that the simu¬ 
lated C-S-H sample was prepared by introducing defects 
inside tobermorite (see section]]). 

For consistent comparison with C-S-H, this system 
was simulated using the REAXFF potential, i.e. the 
same as the one we used for C-S-H, the same time step 
and following the same procedure. We started from an 
initial llA tobermorite cell m composed of 288 atoms, 
relaxed it at zero pressure and 300K during 2.5ns in the 
NPT ensemble, checked the convergence of the volume 
and the energy and run a 25ps simulation in the NVT 
ensemble for statistical average. Note that the use of 
the REAXFF potential does not induce any significant 
change of volume nor structural modification with re¬ 
spect to the starting configuration [16] . 

It should be noted that 11A tobermorite does not have 
the same composition nor the same density as C-S-H. To 
be able to compare the results with C-S-H, we rescaled 
all the computed properties to take this difference into 
account. In practice, this is achieved by replacing the 
concentration of every species by the ones in C-S-H. In 
the following, we will refer to the results obtained with 
the tobermorite system and after the mentioned rescaling 
as being the properties of the crystal. 


Simulation of the glass analogue 

We aimed at comparing the structure of C-S-H with 
the one of an ideal glass. To that end, we created an 
artificial glass analogue of C-S-H by heating and cool¬ 
ing a C-S-H configuration. Note that this methodol¬ 
ogy is commonly used to simulate glassy materials; e.g. 
one usually prepare glassy silica by heating and cooling 
quartz |20| . 

We started from the relaxed configuration of C-S-H. 
The system was then instantly heated at 3000K and re¬ 
laxed during 2.5ns at constant pressure in the NPT en¬ 
semble at this temperature, which is well above the melt¬ 
ing temperature calcium silicate systems (for example 
1813K for the dicalcium silicate crystal). This allowed 
the system to loose the memory of its initial configu¬ 
ration. We checked this by computing the root-mean- 
square displacement of each species and making sure 
that, at the end of the simulation, they were far larger 
than the size of the simulation box. The system was then 
gradually cooled from 3000K down to 300K in the NPT 
ensemble, with a cooling rate q of 20, 40 and 80K/ps. 
During the melting and cooling phases, we imposed a 
pressure of O.IGPa to prevent water molecules to leave 
the system at high temperature, which would lead to ar¬ 
tificially large simulation box. However, this pressure 
is low with respect to the large pressure fluctuations in 
the system due to its small size (typically around lGPa). 
Therefore, we do not expect this pressure to have any sig¬ 
nificant impact on the structure of the simulated glass. 
Once at 300K, the obtained glass was relaxed in the NPT 
ensemble at zero pressure during 2.5ns. At this stage, we 
checked the convergence of the energy and volume of the 
system. Eventually, we run a 25ps simulation in the NVT 
ensemble at 300K for statistical averaging. 

Note that the volume of the simulation box increases 
by 80% while heated at 3000K. However, after being 
cooled, the volume of the system goes back to a value 
fairly close to the volume of the initial C-S-H box (larger 
by 3%). Once again, to allow for a consistent comparison 
with C-S-H, we used the same REAXFF potential and 
the same time step during the entire procedure. In the 
following, we will refer to the results obtained with this 
quenched system as being the properties of the glass. 

SHORT-RANGE ORDER 
Total pair distribution function 

As a first step, we computed the total pair distribu¬ 
tion function (PDF) g^(r) of the three systems from the 
knowledge of the partial PDF gij(r): 

n 

9t ( r ) = c i c j9iji r ) 
i,j =1 


( 1 ) 
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FIG. 1. Total pair distribution function of C-S-H, compared 
with the one of the corresponding glass (top) and of the crystal 
(bottom). 

where Cj is the fraction of i atoms (Si, Ca, H or O). 

Fig. [TJshows the total PDFs of C-S-H, compared with 
the one of the crystal and the glass. We note that the 
PDFs of those three systems does not show significant dif¬ 
ferences that would clearly allow distinguishing a glassy 
from a crystalline phase. They actually present a simi¬ 
lar shape than the PDFs of silicate ED and chalcogenide 
glasses |22| . The first peak in the lA region is associ¬ 
ated to H-0 bonds. The second peak in the 1.6A region 
corresponds to a superposition of the H-H and Si-0 cor¬ 
relations. Following peaks are not easy distinguishable as 
they result of the superposition of several contributions 
of different pairs of atoms. 

Overall, we observe a closer agreement of the PDF of 
C-S-H with the one of the glass than with the one of 
the crystal. This tends to show that C-S-H is mostly 
amorphous. 


Partial pair distribution function 

To better characterize the order and disorder in C- 
S-H, we computed all the partial PDFs. Selected partial 
PDFs are shown in Fig. [2] Note that O w refers to oxygen 
atoms belonging to water molecules. Fig. & and i 3 
show that the short-range order of Si and Ca atoms is 
roughly similar in the crystal, in C-S-H and in the glass. 
In particular, the position and integral of the first peak 
do not show any significant change, which means that 
the atomic bond distance Si-0 and Ca-0 as well as the 
coordination number of the cations are the same. This is 
not surprising as glasses typically tend to retain the same 
local order as observed in the corresponding crystal. 

However, some differences start to appear at larger 
scale. The second coordination shell peak of the Si-0 par¬ 
tial PDF of C-S-H shows a bimodal distribution which 
is reminiscent of the one observed in the crystal. On the 


contrary, the glass only shows a smooth broad peak cor¬ 
responding to the second coordination shell. This means 
that a certain degree of order is maintained in the silicate 
layers in C-S-H. In the case of the Ca-0 partial, the PDF 
of the crystal shows a small peak around 3.2A which dis¬ 
appears in C-S-H. However, the peaks observed in the 
PDF of the glass are broader than the ones of C-S-H. 
Overall, C-S-H seems to present an atomic order that 
is intermediate between the ones of the crystal and the 
glass in the Si and Ca local environment. 

On the contrary, the environment of hydroxyl groups 
(Fig. |2j0 and water molecules show a different be¬ 
havior. Once again, the first peak is fairly similar in C- 
S-H, in the crystal and in he glass, corresponding to sim¬ 
ilar bond distances and coordination numbers. However, 
structural correlations at larger scale (r>1.5 A) appear 
to be almost identical in C-S-H in the glass, and signifi¬ 
cantly differ from the ones observed in the crystal. This 
is a clear evidence the disorder observed inside C-S-H 
mainly arises from water molecules and hydroxyl groups, 
i.e. from the hydration of the network. 

Since the structural differences between C-S-H, the 
glass and the crystal seem to appear at scales larger than 
the one of the first coordination shell, we analyzed sec¬ 
ond neighbor cation-cation and water-water correlations. 
Overall, the first peaks of the Si-Si (Fig. §0 and Ca- 
Ca (Fig. §>) partial PDFs of C-S-H are systematically 
broader than in the crystal, but sharper than in the glass. 
This confirms the fact that the calcium and silicate lay¬ 
ers in C-S-H retain some crystal-like order, as opposed 
to H-H (Fig. &) and water-water (Fig. i*) correlations 
which are essentially amorphous. 


MEDIUM-RANGE ORDER 

As we realized that the local environment of the dif¬ 
ferent species was not enough to distinguish the order in 
C-S-H from the one of the crystal and the glass, we an¬ 
alyzed the medium range order (MRO) of those systems 
by computing the neutron structure factor. The partial 
structure factors were first calculated from the pair dis¬ 
tribution functions g %3 (r): 

S ij (Q) = l+Q 0 f 4tt r 2 {gij (r)-1) sm ^ r ) p L (r) d?- (2) 
Jo 

where Q is the scattering wave vector, g 0 is the aver¬ 
age atom number density and R is the maximum value 
of the integration in real space (here R = 6A). The 
Fi,(r) = sin(7rr/ R )/ (nr / R) term is a Lortch-type win¬ 
dow function used to reduce the effect of the finite cutoff 
of r in the integration [22. As discussed in the reference 
[2T , the use of this function reduces the ripples at low Q 
but induces a broadening of the structure factor peaks. 
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FIG. 2. Partial distribution functions for the Si-0 (a), Ca-0 (b), H-O (c), and O w -H pairs, where O w refers to oxygen atoms 
in water molecules. For each pair, the partial distribution function of C-S-H is compared with the one of the corresponding 
glass (top) and of the crystal (bottom). 



FIG. 3. Partial distribution functions for the Si-Si (a), Ca-Ca (b), H-H (c), and O w -O w pairs, where O w refers to oxygen atoms 
in water molecules. For each pair, the partial distribution function of C-S-H is compared with the one of the corresponding 
glass (top) and of the crystal (bottom). 
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FIG. 4. Total neutron structure factor of C-S-H, compared 
with the one of the corresponding glass (top) and of the crystal 
(bottom). 


The total neutron structure factor can then be evaluated 
from the partial structure factors following: 


^ ^ CiCjbibj) ^ ' CiCjbibjSij(Q) (3) 
i,j =i i>i = i 

where Cj is the fraction of i atoms (Si, Ca, H or O) and 6,; 
is the neutron scattering length of the species (given by 
4.1491, 4.70, -3.7390 and 5.803fm for silicon, calcium, hy¬ 
drogen and oxygen atoms respectively [25]). Once again, 
note that corrections are applied to take into account the 
different composition and density of C-S-H with respect 
to tobermorite. 


Neutron structure factor 

Kg. a shows the total neutron structure factor for C- 
S-H as well as for the corresponding crystal and glass. 
Once again, the three structure factor do not show sig¬ 
nificant differences and are typical of the ones observed 
in silicate [2B] and chalcogenide m glasses. It should be 
mentioned that the limited sizes of the simulated systems 
do not allow us to study large-scale correlations (large r, 
low Q). Focusing on the low Q region of the structure 
factor, we note that the peaks show the same position in 
C-S-H and in the glass, although they appear sharper in 
C-S-H, signature of an increased order on the MRO. On 
the contrary, the peaks at low Q do not show the same 
positions. Since each peak is a signature of a typical spa¬ 
tial repetition distance in the MRO m, it appears that 
the atomic order at intermediate length scale of C-S-H is 
amorphous. In particular, we observe a sharp first peak 
around 1.5A -1 for the crystal, resulting of a strong order 
at large r, which is absent in C-S-H and the glass. 


Partial structure factors 

Similarly to the case of the total PDF, the apparent 
disorder observed in the total structure factor actually 
results of some order and some disorder for different pair 
of species. This is what we investigated by studying the 
decomposition of the total structure factor into the con¬ 
tributions of each partial (see Fig. [5] and [g]) . Conclusions 
actually appear to be the same of the ones of section [] 
The MRO observed in calcium and silicate layers (tracked 
in the Si-O, Ca-O, Si-Si and Ca-Ca partials) for C-S-H 
is closer to the one of the crystal. In particular, Si-O, 
Ca-0 and Si-Si partials respectively show peaks, respec¬ 
tively around 3.5, 4.5and 3.5 A -1 , for both C-S-H and 
the crystal, which are absent in the case of the glass. On 
the contrary, partial associated to water and hydroxyl 
groups are almost indentical for C-S-H and the glass, 
and significantly differ from those observed in the crystal. 
Hence, the disorder observed in the medium range order 
of C-S-H is mainly driven by those hydration species. 


CONCLUSION 

A detailed comparison of the structure of C-S-H as 
well as the one of corresponding crystal and glass has al¬ 
lowed us to better characterize the atomic order inside 
the binding phase of cement. Overall, it appears that 
the structure of C-S-H is closer to the one of a glass 
than to the one of a crystal. However, some atomic or¬ 
der reminiscent from the one of a crystal is still found. 
This manifests by a fairly layered structure, by silicate 
chains that are not completely amorphous and by a non- 
random distribution of calcium atoms. On the contrary, 
water molecules and hydroxyl groups show a completely 
glassy spatial distribution, which suggests that the over¬ 
all disorder is mainly caused by its hydration. 

Hence, the atomic order of C-S-H can be qualified 
as being intermediate, definitely not fully crystalline but 
not perfectly amorphous. The combination of a layered 
structure with an overall disorder may explain the para¬ 
doxical observation of an ordered or disordered structure 
according to the experimental technique that is used m- 
Indeed, when averaged on all pairs of atom, C-S-H ap¬ 
pears to be fully glassy. Thus, only the experiments pro¬ 
viding an access to the detail of the order and disorder 
around each species would be able to capture the com¬ 
plexity of C-S-H. 


* Contact: bauchy@ucla.edu; 

Homepage: littp://mathieu.bauchy.com 
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FIG. 5. Partial structure factors for the Si-0 (a), Ca-O (b), H-O (c), and O w -H pairs, where O w refers to oxygen atoms in 
water molecules. For each pair, the partial structure factor of C-S-H is compared with the one of the corresponding glass (top) 
and of the crystal (bottom). 
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FIG. 6. Partial structure factors for the Si-Si (a), Ca-Ca (b), H-H (c), and 0 W -0 W pairs, where 0 W refers to oxygen atoms in 
water molecules. For each pair, the partial structure factor of C-S-H is compared with the one of the corresponding glass (top) 
and of the crystal (bottom). 
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